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1. Introduction 

This paper is a continuation of three different stories. The first story begins with the work 
of Erdos and Renyi [19] on the random graphs Q{n,p) (with n vertices, and where every 
edge is included independently with probability p ), and their connectivity properties (see 
also [12]). It continues with Penrose’s work [40] extending this to the geometric graphs 
Q{V,r) - graphs whose vertices are random point processes {V) in a d-dimensional space, 
and whose edges are determined by the proximity between the points (distance less than 
r). In both models, a phase transition occurs when the number of vertices n goes to 
infinity, and the average degree is of order logn {p = or r oc In this 

sharp phase transition, if the degree is lower then a known threshold the graph has many 
connected components, but once the threshold is passed the graphs become connected. 

The second story begins much later in the work of Linial and Meshulam [31], and is 
about high dimensional extensions of the Erdos -Renyi theory, where instead of a graph 
one studies a random simplicial complex. See [4, 5, 6, 17, 18, 26, 28, 29, 31, 32] for a 
variety of important results. Random simplicial complexes have various models, and much 
is still unknown about those, although remarkable phenomena have been found regarding 
their homology (a higher dimensional generalization of connectedness, see Section 2.1) 
and their fundamental groups. A direct motivation for the work presented in this paper 
can be found in Kahle’s work [26] describing the vanishing of homology of random flag 
complexes. Random flag complexes are generated from Erdos -Renyi graphs by adding 
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a fc-simplex for any {k + l)-cliques in the graph. The main result in [26] states that the 
fc-th homology group vanishes (becomes trivial) when p = ^ ^ and the 


result of Erdos and Renyi can be viewed as the special case fc = 0. Our goal in this paper 
is to explore an equivalent extension in the random geometric graph context. 

The third story is related to the theory of coverage processes (cf. [23]). In particular, 
Flatto and Newman [20] have shown that on a d-dimensional unit volume Riemannian 
manifold A4, generating n uniformly distributed points and covering them with balls of 
radius r will cover M with high probability provided that r oc (iHii^i)i/'^ (resembling 
the connectivity threshold discussed above). On first sight, this phase transition might 
seem unrelated to the problem at hand, but as we shall see soon, it provides significant 
information about the vanishing of homology in geometric complexes. 

There are various ways to dehne geometric simplicial complexes. The one studied in 
this paper is the random Cech complex. In this model we generate points randomly in a 
metric space (e.g. a compact Riemannian manifold A4), and form a simplicial complex 
by hxing a radius r and asserting that k points span a (fc — l)-simplex if the r-balls 
around these points have a nonempty intersection. The theory of random geometric 
complexes (such as the Cech complex) has been growing rapidly in the past decade (cf. 
[1, 2, 9, 11, 25, 27, 39, ?, 41] and the comprehensive survey in [10]). In particular, the 
limiting behavior of the homology of these complexes has been studied, when n —>■ oo 
and r = r(n) —>■ 0. 

In this paper we study the following model. Let be the d-dimensional 

flat cubical torus (see Section 3), and let Vn be a homogeneous Poisson process in 
with intensity n (i.e. E{|7^„|} = n). We are interested in phase transitions related to 
the homology of the Cech complex C{'Pn,r) as n —>■ oo and r = r(n) —)■ 0. The limiting 
behavior of the complex is controlled by the term A := w^nr'^, where uJd is the volume 
of a d-dimensional Euclidean unit ball. This quantity can be viewed as measuring the 
expected degree of the underlying graph. 

For this random geometric complex, we shall see that essentially there are two sharp 
phase transitions. The first transition is the Erdos-Renyi-Penrose connectedness type, 
which occurs when A = 2“‘^logn. The second transition is related to the rest of the 
homology groups, and occurs when A = logn. Below the second threshold, we see a 
lot of extra homology (cycles, or holes in various dimensions k < d), while above the 
threshold the homology of the random complex becomes identical to that of the torus (a 
phenomenon we refer to as the “vanishing” of homology). 

As mentioned earlier, a phenomenon related to the vanishing phase transition has also 
been observed in the literature on coverage processes (see [23]). Flatto and Newman [20] 
have shown that A = logn -f (d — 1) log logn is the threshold for the balls or radius r 
around Vn to cover a d-dimensional manifold (see also Appendix B). Since the union 
of balls and the Cech complex share the same homology (see Section 2.2), and since 
coverage implies the vanishing of homology, the result in [20] provides an upper bound 
to the vanishing threshold we are interested at. 

Our main results (see Section 5) assert that when A = (1 -f e) logn the homology of 
C{Vn,r) is identical to that of while prior to this point, i.e. when A = (1 — e) logn. 
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the homology groups of C{Vn,r) are very large (the number of cycles grows to infinity). 
To prove this statement we show that at the moment before coverage many nontrivial 
cycles are still being formed. We do this using Morse theory for distance functions (see 
e.g. [16, 21], and Section 2.3). In Morse theory (see Section 2.3.2), one counts critical 
points of functions according to the index of the Hessian (i.e. the matrix of second order 
derivatives) at the critical points, with an aim towards computing homology. The main 
idea is that a critical point of index k may either give rise to a /c-dimensional cycle or 
alternatively signal the death of a (fc — l)-dimensional cycle. In general, however, there 
is no rule to determine which of these takes place. As a result, Morse theory ordinarily 
does not provide an accurate analysis of homology, but rather a set of inequalities. In this 
paper on the other hand, we construct a special type of critical points that are guaranteed 
to generate cycles in the dimension of their index. We refer to these cycles as “0-cycles” 
(see Section 7), and use them to show the existence of non-trivial homology. 

The study of critical points leads to some additional information about the finer scale 
structure. We see that analogously to the results of [20], at an additive log log n scale, 
there seems to be a separation between the vanishing of the different dimensions of 
homology. In particular, our results (see Theorem 5.4) suggest that the exact vanishing 
threshold of the homology in dimension k is in the range 

[log n + {k — 2) log log n, log n + k log log n]. 

For example, at this fine scale the 1-dimensional cycles vanish before the 3-dimensional 
ones (with high probability). It is very interesting but beyond the scope of our current 
methods to try to understand exactly where in these finer scale ranges the various ho¬ 
mology groups converge to those of the torus, and the nature of the vanishing results 
that might occur at that scale. 

While random geometric complexes provide grounds for rich and deep theoretical prob¬ 
ability research, we note that they also find applications in data analysis and network 
modeling. The rapidly developing field known as Topological Data Analysis (TDA) fo¬ 
cuses on using topological signatures of data in machine learning and statistics (for some 
introduction see [14, 22, 42]). Geometric complexes play a key role in the conversion of ab¬ 
stract topological questions into a simplified set of algebraic and combinatoric operations 
that can be coded in software. Analyzing the behavior of random geometric complexes 
is therefore imperative in order to provide TDA with rigorous statistical statements (see 
e.g. [15, 37, 38]). The results in this paper provide information that is asymptotically 
an improvement of the few estimates on the sample complexity in topological inference 
problems. It is also related to [7, 8] on the rate of convergence for homology estimation 
in riemannian manifolds. 

The structure of this paper is as follows. Sections 2-4 provide the necessary background 
and definitions for this work. The main results are presented in Section 5, describing the 
phase transition for the vanishing of homology. Sections 6 and 7 present the main ideas 
and lemmas used to prove the upper and lower bounds, respectively. Section 8 presents 
the detailed proofs for the statements in Sections 5-7. 

Finally, while the results of this paper are proved explicitly only for the case of the 
cubical torus, they trivially apply to any flat torus (of unit volume, i.e. K'^/L for any 
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lattice L in In section 9 we give a heuristic (ineffective) argument why they also 
should apply to arbitrary closed Riemannian manifolds (normalized to have unit volume). 
We hope to return to this question with effective estimates in a later paper. 


2. Topological background 

In this section we wish to provide a brief introduction to homology, Cech complexes and 
Morse theory, which we will use later in this paper. 


2.1. Homology 


We wish to introduce the concept of homology here in an intuitive rather than a rigorous 
way. A comprehensive introduction to the topic can be found in [24, 36]. Let X be a 
topological space, the homology of X is a sequence of abelian groups denoted {lIi(X)}“g. 
Homology is a topological invariant, namely if / : X —>■ X is a homeomorphism, then it 
induces an isomorphism /* : H^{X) H^{Y) between the homology groups. 

In the case where homology is computed using coefficients in a held F, then Hi{X) is 
simply a vector space^. The basis elements of zeroth homology Hq{X) correspond to the 
connected components of X. For example, if X has three connected components, then 
Hq{X) = (where = denotes isomorphism), and each of the basis element corresponds 
to a different connected component of X. For fc > 1, the basis elements of the fc-th 
homology Hk{X) correspond to fc-dimensional “holes” or (nontrivial) “cycles” in X. An 
intuitive way to think about a fc-dimensional cycle is as the boundary of a (fc -I- 1)- 
dimensional body. For example, if X a circle then Hi{X) = F, if X is a 2-dimensional 
sphere then H 2 {X) = F, and in general if X is a d-dimensional sphere, then 


Hk{X) 


F fc = 0, d 
{0} otherwise. 


Another example which will be relevant in this paper is the torus. The 2-dimensional 
torus (see Figure 3) has a single connected component and a single 2-dimensional 
hole (the void inside the surface). This implies that ido(T^) = F, and = F. As 

for I-cycles (or closed loops) the torus has two different loops, and so idi(T^) = F^. In 
this paper we will consider the general d-dimensional torus denoted by In that case 
it turns out that idfc(T‘^) = f(*=). 

The last term we want to introduce here is homotopy equivalence. Let X, Y be two 
topological spaces. A homotopy is a continuous function F : X x [0,1] —)■ X, which can be 
viewed as a indexed sequence of functions F{-,t) : X Y. Two functions fo, fi : X —>■ X 
are called homotopic (denoted /q ~ /i) if there exists a homotopy F such that F{-, 0) = fo 
and F(-,l) = fi. Finally, / : X —>■ X is called a homotopy equivalence if there exists 

^We introduce homology with field coefficients for simplicity. Our results, however, apply to homology 
with arbitrary coefficients, not just fields. 
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g : Y ^ X such that g o f = idx and / o g ~ idy (idx refers to the identity function 
on X). If there exists a homotopy equivalence f : X ^ Y then we say that X and Y are 
homotopy equivalent, and denote it by X ~ y. Loosely speaking, X ~ y means that we 
can continuously transform one of the spaces into the other, generalizing the notion of 
homeomorphic spaces (for example, a ball and a point are homotopy equivalent but not 
homeomorphic). In particular, if X ~ y then it can be shown that they have the same 
homology, i.e. iLfe(X) = Hk{Y) for all k. A space that is homotopy equivalent to a single 
point is called contractible. 


2.2. Cech complexes 

For a given set S, an abstract simplicial complex S on S' is a collection of finite subsets 
A S 2'® such that if A € S and B <Z A then also B G Y. We refer to sets A G E 
with |A| = fc + 1 as the fc-simplexes or fc-faces in E. In this paper we study an abstract 
simplicial complex known as the Cech complex, defined in the following way. 

Definition 2.1 (Cech complex). Let V = {xi,X 2 , ■ ■ ■ ,Xn} be a collection of points in 
a metric space (X,p), and let r > 0 and let Br{x) be the ball of radius r around x. The 
Cech complex C{V,r) is constructed as follows: 

1. The 0-simplexes (vertices) are the points in V. 

2. A k-simplex [xi^,..., Xi^ ] is inC(V,r) if [fj^^Brixi.) ^ 0. 



Figure 1. A Cech complex generated by a set of points in The complex has 6 vertices (0-simplexes), 
7 edges (1-simplexes) and one triangle (a 2-simplex). 

Associated with the Cech complex C(P,r) is the union of balls used to generate it, 
which we define as 

U{V,r)=[j Br{p). (2.1) 

p&V 

The spaces C{V,r) and U{V,r) are of a completely different nature. Nevertheless, the 
following lemma claims that topologically they are very similar. This lemma a special 
case of a more general topological statement originated in [13] and commonly referred to 
as the ‘Nerve Lemma’. 
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Lemma 2.2. Let C{V,r) andU{V,r) as defined above. If for every Xi^, ... ,Xif. the 
intersection Br(xi^) n • • • n Br{xi^) is either empty or contractible (homotopy equivalent 
to a point), then C{V,r) ~U{V,r), and in particular, 

Hu{C{V,r))^Hk{U{V,r)), Vfc > 0. 

Consequently, we will sometimes be using U{fP,r) to make statements about C{V,r). 
This will be very useful especially when coverage arguments are available. Note that in 
Figure 1 indeed both C(7^, r) a.ndU{V, r) have a single connected component and a single 
hole. 

2.3. Morse theory for the distance function 

The distance function defined below and its critical points will play a key role in our 
analysis of the vanishing threshold for homology. In this section we wish to define the 
distance function and its critical points, and briefly introduce Morse theory. 

Let 7^ be a finite set of points in a metric space {X, p). We define the distance function 
from V as follows - 

Pj,{x) = mmp{x,p), Vx G X. (2.2) 

pev 

Our interest in this function stems from the following straightforward observation about 
the sublevel sets of the distance function: 

P(:^((-oo,r]) = p~^{[0,r]) = U{V,r). 

In other words, the sublevel sets of the distance functions are exactly the union of balls 
associated with the Cech complex. The idea of Morse theory is to link the study of critical 
points of functions with the changes to the homology of their sublevel sets. In analyzing 
the vanishing of homology for the Cech complex, this link will show up to be highly 
useful. 

2.3.1. Critical points of the distance function 

The classical definition of critical points in calculus is as follows. Let / : —>■ K be a 

function. A point c S M is called a critical point of / if X f{c) = 0, and the real number 
/(c) is called a critical value of /. A critical point c is called non-degenerate if the Hessian 
matrix Hf(c) is non-singular. In that case, the Morse index of / at c, denoted by p{c) is 
the number of negative eigenvalues of Hf{c). 

Note that the distance function defined in (2.2) is not everywhere differentiable, 
therefore the definition above does not apply. To overcome this problem we will use Morse 
theory for min-type functions which was developed in [21]. Min-type functions are of the 
form: f{x) = mini<i<m Q.iix) where oi,..., am are smooth functions. The work in [21] 
provides the definitions of regular and critical points for this type of functions, and shows 
that large parts of Morse theory (discussed below) can be applied to these functions as 
well. 
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In this section we wish to introduce the definitions of critical points for for the 
special case where V (with the Euclidean metric p{x,y) = \\x — j/||). Later on, we 

will see how to extend these ideas to the d-dimensional torus. Note that the distance 
function p^ does not fall directly into this category of min-type functions, since the 
function p(p, •) is non-differentiable at p. However, the squared-distance function is 
a true min-type function for which the results in [21] apply. The sublevel sets of are 
still union of balls, i.e. (Pp)“^((—c»,rj) = U{'P,yC-), and therefore any conclusion we 
can make using Morse theory for p^ could be naturally translated to similar statements 
about the sublevel sets of p^. Hence, from here on, we will talk about critical points and 
Morse theory for p^, while formally we mean p^. 

The following definitions for the critical points of p^ appeared already in [9], and 
we include them here for completeness. They are merely an adaptation of the general 
definitions appeared in [21] (Section 3.1) to the distance function case. Start with the 
local (and global) minima of p^, which are all the points in V where p^ = 0, and call 
these critical points with index 0. For higher indices, we have the following. 

Definition 2.3. A point c is a. critical point of p^ with index 1 < k < d if there 
exists a subset y <ZV of k -\- 1 points such that: 

1. Vy ey : Pt,{c) = \\c-y\\, and, Vp G V\y we have \\c-p\\ > Pp(p). 

2. The points in y are in general position (i.e. the k + 1 points of y do not lie in a 
(fc — 1)-dimensional affine space). 

3. The point c lies inside the open k-simplex spanned by y. 

The first condition ensures that p^ = p^ in a neighborhood of c. The second condition 
implies that the set y lies on a unique (k — 1)- dimensional sphere. We note that this 
generality requirement is satisfied almost surely when the points are generated by a 


Poisson process. We shall use the following notation: 

5(3^) = The unique (fc — l)-dimensional sphere containing y, (2.3) 

C{y) = The center of S{y) in K'', (2.4) 

i?(3;) = The radius of 5'(3^), (2.5) 

B{y) = The open ball in with radius R{y) centered at C{y), and (2.6) 

A(3^) = The open A:-simplex spanned by the k-\-\ points in y. (2.7) 


Note that 5'(3^) is a (fc — l)-dimensional sphere, whereas B{y) is a d-dimensional ball. 
Obviously, S{y) C dB{y), but unless k = d, S{y) is not the boundary of B{y). Since 
the critical point c in Definition 2.3 is equidistant from all the points in y, we have that 
c = C{y). Thus, we say that c is the unique index k critical point generated by the k-\-l 
points in y. The last statement can be rephrased as follows: 

Lemma 2.4. A subset y GV of k -\-1 points in general position generates an index k 
critical point if and only if the following two conditions hold: 

(1) C(3^)gA(37), and (2) P n 5(3^) = 0 
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The critical point generated is C{y), and the critical value is R{y) = Pj,{C{y)). 

Figure 2 depicts the generation of an index 2 critical point in by subsets of 3 points. 

Vs 



Figure 2. Generating a critical point of index 2 in i.e. a maximum point. The small blue disks are the 
points of V. We examine three subsets of V: ^2 = ^3 = 

S{yi) the dashed circles, whose centers are C{yi) = Ci. The shaded balls are B(yi), and the interior 
of the triangles are A(y’i)' The arrows represent the flow direction. For both conditions in Lemma 
2.4 hold and therefore ci is a critical point. However, for 3^2 condition (1) does not hold, and for 3^3 
condition ( 2 ) fails, therefore C 2 and C 3 are not critical points. 


2.3.2. Morse Theory 

The study of homology is strongly connected to the study of critical points of real valued 
functions. The link between them is called Morse theory, and we shall describe it here 
briefly. For a more detailed introduction, we refer the reader to [34]. 

The main idea of Morse theory is as follows. Suppose that is a closed manifold (a 
compact manifold without boundary), and let / : A4 —>■ K be a Morse function. Denote 

Mv := f~^{{—oo,v]) = {x G M : f{x) < v} C M 

(sublevel sets of /). If there are no critical values in (a, b], then Ma and Mb are homotopy 
equivalent and in particular Hk{Ma) — Hk{Mb) for all k. Next, suppose that c is a 
critical point of / with Morse index k, and let v = /(c) be the critical value at c. Then 
at Mv homology changes in the following way. For a small enough e we have that the 
homology of Mv+e is obtained from the homology of Mv-e by either adding a generator 
to iJfc (increasing its dimension by one) or removing a generator of (decreasing its 

dimension by one). In other words, as we pass a critical value, either a new fc-dimensional 
cycle is formed, or an existing (k — l)-dimensional cycle is terminated. 

While classical Morse theory deals with smooth (or C^) Morse functions on compact 
manifolds [34] it has many generalizations, and the extension to “min-type” functions 
presented in [21] enables one to apply similar concepts to the distance function p^, as 
well. In particular, the critical points defined in Section 2.3.1 have similar effect on the 
homology as in classical Morse theory. 
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The results in this paper should apply to homogeneous Poisson processes generated on 
any compact Riemannian manifold, as we discuss in Section 9. However, for the sake 
of this first investigation we will focus on a single special case - the d-dimensional flat 
cubical torus. 

By ‘flat torus’ we refer to the quotient In other words, the flat torus can 

be thought of as the unit cube in with its opposite sides “glued” together (see Figure 
3). Consequently, the metric p on T'’*, known as the toroidal metric, is given by 

p{x, y) = min \\x - y + A|| , x,y & [0,1]''. 


The main advantage of the flat torus is that it allows us to work with a simple (almost) 
Euclidean metric, while avoiding boundary effects that exist, for example, in the cube. 



"It- 


Figure 3. The flat torus is obtained by taking the unit cube in and identifying the two pairs of 
opposite sides. This identification creates the two loops in the torus. 


More specifically, the radius of convexity rconv of a Riemannian manifold M is dehned 
as the largest value r such that every ball of radius r in is convex. It is also the 
threshold below which any two points in a Riemannian ball can be connected by a 
unique short geodesic that lies entirely within the ball. In the case of flat torus we have 
fconv = 1/2, and it can be shown that every ball in of radius less than r^onv can be 
isometrically embedded into Since we will mostly be considering infinitesimally small 
neighborhoods, our calculations will be as if we are working in the Euclidean space. 

Our choice to state and prove the results in this paper for the flat torus (T'^, p) stems 
from the relative elegant calculations involved using the toroidal metric. While locally, the 
torus (as any other closed manifold) looks almost identical to the Euclidean space 
some caution is required when considering global phenomena. Next, we wish to discuss 
the necessary adjustments needed to study the Cech complex and the distance function 
presented in Sections 2.2 and 2.3. 


3.1. The Cech complex on the torus 

A large part of our analysis will rely on the Nerve Lemma 2.2 (stating that C{'P,r) ~ 
U{V,r)). For this lemma to apply to C{V,r) we need the intersections of balls of radius 
r to be either empty or contractible. For balls in this is always true, but on the torus 
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(and other compact manifolds), this is not true for all radii. For example, on the torus, 
even a single ball of radius 1/2 is not contractible (it covers the 1-cycles of the torus). It is 
true however, for any radius smaller than rconv, and therefore we limit our discussion to 
the Cech complex C(P,r) with r < rconv As we will see later, the vanishing of homology 
which we are seeking occurs with high probability at r = o(rconv)i and therefore this 
restriction will not affect our results. 


3.2. The distance function on the torus 

A similar caution is required studying the Morse theory for the distance function. In 
Section 2.3 we defined the notion of critical points for the distance function for 
7^ C using the Morse theory for min-type functions in [21]. We were able to do so 
due to the fact that p^(p, •) : —>■ K is smooth, and therefore = min^gp p^(p, •) is 
a min-type function. However, on the torus (as any compact Riemannian manifold) the 
situation is more delicate. If p G then /5^(p, •) : —>■ M is smooth at x G only 
if the geodesic from x to p is unique. Consequently, we have to be careful about how 
we apply Morse theory for min-type function in this case. We start with the following 
lemma. 

Lemma 3.1. Let V C T'^ be a finite set. If r < rconv then : U{'P,r) R is a 
min-type function. 

Proof. Although p^(p, •) is not everywhere smooth in it is smooth in the ball 
Br{p) for every r < Xconv Let x G be such that Pp(x) < rconv Denoting '■= 
{per : p(p,x) < ^conv} we can find a small neighborhood such that p(p, p) rconv 
for every p G Vx and y G Nx and such that 

Pv iv) = min p(p, y), Vy G Nx. 

pe'Px 

Thus, in Nx we have that is a min-type function (the minimum of smooth functions). 
We conclude that if r < Xconv then in the set U{V,r) the function p^ defines a germ of 
a min-type function. 

□ 

Remark: This lemma allows us to study the homology of U{V^ r) for r < Xconv using 
the distance function, similar to Note that the restriction r < rconv is similar to the 
one that arose form the Nerve Lemma. 

Another issue we need to resolve is the uniqueness of critical points. If 7^ C we saw 
in Lemma 2.4 that the critical point generated by a subset 3^ C 7^ is unique. On the torus, 
this is not necessarily true, even if r < Xconv For example, the set V = {(0.2,0), (0.8,0)} C 
generates two critical points of index 1 (saddle points) at ci = (0,0) and C 2 = (0.5,0), 
with Pp(ci) = 0.2 and Pp(c 2 ) = 0.3, both values are less than rconv = 0.5. Requiring that 
r < rconv/3 resolves this issue as the following lemma states. 
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Lemma 3.2. Let r < rconv/3. If c G U{V,r) is a critical point of index k, then it is 
generated hy a subset y G V of k + 1 points (similarly to Definition 2.3). Furthermore, 
c is the unique critical point inU{V,r) generated by y. 

Proof. For every x G UifP, r) note that the values of in the ball Br{x) are completely 
determined by points in the set V n Bsrix), i.e. Isr-ix) = \b,.{x)- If r < rconv/3 

the ball B^r{x) can embedded isometrically in which implies that p^|_B^(a;) behaves 
identically to the distance function in discussed in Section 2.3. If c G U{V,r) is a 
critical point of p^, then using the isometry i? 3 r(c) —)■ Definition 2.3 applies to it 
as well. In particular, following Lemma 2.4, c is generated by a subset 3^ C 7^ of fc + 1 
points, such that Pj,(c) = Py{c) < r. 

Suppose that there exists another critical point c! G U{V,r) generated by y. Since 
Pj,{c') = Py{c') < r as well, we must have d G i?2r(c) C B^j.{c). However, since in 
the set y generates a unique critical point, the same should hold for the embedding 
Bsrid) —>■ and therefore c = d. 

□ 

Remark: Using the uniqueness of critical points in Ui(P,r) (r < rconv/3), we can 
define S(y), C(y), R(y), B{y), and A(3^) similarly to (2.3)-(2.7) using the local isometry 
described above. 

To conclude, while the distance function on the torus T”* introduces some features 
that do not occur in (e.g. singular points, and non-uniqueness of critical points), if 
r < rconv/3 then in U{V,r) C the critical points of the distance function behave the 
same was as in Consequently, we define 

rmax ■— rconv/3, (3-1) 

which will be useful for us later. As mention earlier, the vanishing threshold we look for is 
much smaller than rconv (and rconv/3). Therefore, we will be able to analyze the distance 
function on the torus similarly to the Euclidean space. 


4. Definitions and notation 

4.1. The Poisson process 

In this paper we study the homology of random Cech complex constructed from a random 
set of points. The points are generated by a homogeneous Poisson process which can be 
defined as follows. Let (X,p) be a compact metric-measure space and let Ai, A 2 ,... be 
a sequence of i.i.d. random variables uniformly distributed in X. Let N ^ Poisson (n) be 
a Poisson random variable, independent of the A^-s. Then we define 

iP„ = {Ai,...,A^}. (4.1) 

This definition is equivalent to saying that is a homogeneous Poisson process with 
intensity n. In particular then the following holds: 
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1. For every compact measurable set A: {Vn n A| ^ Poisson {nfi{A)) 

2. If A and B are disjoint, then the variables Vn<^ A and Vn H B are independent 
(this property is usually referred to as ‘spatial independence’). 

We will be interested in asymptotic behavior as n —>■ oo. For any event A that depends 
on Vn, we say that A occurs with high probability (w.h.p.) if lim„_>oo P (^) = 1- 

4.2. Some notation 

The following objects will be used repeatedly in the paper: 

• The union of balls - U{n, r) :=U{Vn,r), 

• The Cech complex - C(n,r) := C{Vn,r), 

• The distance function - p„(-) := 

Throughout the paper, k will mostly be a fixed positive integer value referring to 
either an index of a critical point or degree of homology. For a fixed k, we will use y 
to represent a set of fc + 1 random variables. When evaluating probabilities or moments 
using integrals will use the following notation - 

• a: G T'^ - a single variable, 

• X = (xo,..., Xk) G - a (A: + l)-tuple of points in 

• y = (yi, ■ ■ ■, Vk) G - a fc-tuple of points in M”*, 

• If / : (Td.)>^+^ M then /(x) := f{xo,Xk), 

• If / : K then /(O, y) := /(O, yi,...,yk). 

In random Cech complexes, similarly to random geometric graphs, the term that 
controls much of the limiting behavior of the complex is 

A := ujdnr^, (4.2) 

where cod is the volume of a d-dimensional Euclidean unit ball. The value A is the expected 
number of points in a ball of radius r. This value is directly related to the expected vertex 
degree (number of neighbors) in the random Cech complex C(n, r). The higher A is, the 
denser the complex is. The vanishing thresholds we study in this papers are values of A 
that guarantee the vanishing of homology. 


Finally, we will 

use 

the following asymptotic notation. 

• a(n) s 

s b{n) 


lim„_,.oo ain)/b{n) = I, 

• a(n) ^ 

- bin) 


lim„^oo a(n)/&(n) G (0,oo), (i.e. a(n) = 6(&(n))) 

• a(n) <C b{n) 


lim„_>oo a{n)/b{n) = 0 (i.e. a(n) = o(&(n))). 


5. Main results 

Let Vn be the homogeneous Poisson process in with intensity n, and C{n,r) be the 
corresponding Cech complex. If r is large enough we expect the union of balls U (n, r) to 
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cover the torus and then, by the Nerve Lemma 2.2 we have Hk{C{n,r)) = for 

every 0 < k < d. Translating the results in [20] to our case (see Corollary B.2) we have 
that the threshold value for coverage is 


Ac := log n+ {d — 1) log log n. 


In addition, the analysis of connectivity in random geometric graphs (cf. [40]) yields that 
the threshold for connectivity of C{n,r) is 


Ao 


logn 

2 ^ 


(5.1) 


In other words, Ag is the threshold value above which we have HQ{C{n,r)) = 
Analogously to other models of random simplicial complexes (e.g. [26, 33]), we expect 
that between Ag and Ac there would be an increasing sequence of threshold values 


Ag < Ai < Aa < ... < Ad < Ac, 

such that Afc is the threshold value to have Hk{C{n, r)) = We refer to Ak as the 

‘vanishing threshold’ for Hk, in the sense that above this value we see only the fc-cycles 
that belong the torus, and all the other cycles that appeared in lower radii vanish. 

We need to separate the top (d-dimensional) homology from the rest of the homology 
groups. Prior to coverage, since U{n,r) is a “nice” (compact and locally contractible) 
proper subset of T'’* then necessarily Hd{iU{n,r)) = 0 (cf. [24]). Once we cover the torus, 
we know that U{n,r) = T'^, and therefore Hd{U{n,r)) = Corollary B.2 then 

immediately implies the following. 


Theorem 5.1. Suppose that w{n) —)■ oo as n ^ oo. 


lim F{Hd{C{n,r))^Hd{T‘^)) 


1 A = logn-I-(d — 1) loglogn-I-w(n) 
0 A = logn-|-(d — 1) loglogn — w(n). 


In other words, A^ « logn, while Ag « 2“'’* logn. Comparing Ad and Ag (5.1) a 
reasonable conjecture would be that the other vanishing thresholds are of the form A^ « 
Ck logn for some Ck S (2““*, 1). The main theorems of this paper provide bounds for Ak 
{k = 1,... ,d — 1) and show that this conjecture is actually wrong and in fact A^ « logn 
for all 1 < fc < d. 

To analyze the vanishing of the k-th homology (l<fc<d—1), we start by bounding 
the expected value of /3k{r) - the fc-th Betti number of C(n,r) (i.e. the number of non¬ 
trivial fc-cycles). 


Proposition 5.2. // r —?► 0 and A oo, then for every 1 < k < d — 1 there exist 

constants Ok ,bk > 0 such that 

OkuA^-^e-^ < E{/3k{r)} < /3fc(T‘^) + bkuA^e-^. 
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Consequently, we have the following phase transition. 


Corollary 5.3. Suppose that w{n) —)■ oo as n —>■ oo and let 1 < k < d — 1. Then 


lim E{&M} = A = log„ + fcloglog.. + „.(„) 


OO A = logn + (fc — 2) loglogn — 


Finally, we would like to turn the statement about the expectation into a statement 
about the vanishing probability. This results in the following theorem. 


Theorem 5.4. Suppose that w{n) —)■ oo as n ^ oo and let 1 < k < d — 1. 



A = log n + k log log n + w{n) 

A = logn + {k — 2) log log n — w{n), 


where for the lower threshold we require further that w(n) ^ log log logn. 


Theorem 5.4 suggests that the vanishing threshold for the k-th homology satisfies 


Afc G [log n+ {k — 2) log log n, log n + fc log log n]. 
In particular we have the following straightforward corollary. 


Corollary 5.5. For any 1 < k < d — 1, and for any e G (0,1), 



A = (1 + e) logn 
A = (1 — e) logn. 


This indicates that the conjecture made above about the vanishing thresholds is wrong. 
Except for Hq, all other homology groups Hk vanish at the same scale A « logn (i.e. all 
the Cfc, fc > 1 coincide), and the differences in the threshold values show up as a second 
order log logn term. This result is highly counterintuitive, and we consider it one of the 
main contributions of this paper. In Section 10 we will provide an intuitive explanation 
for this phenomenon. 

We provide the proofs for the theorems above in Section 8. The arguments required to 
prove the upper and lower bounds are different, but both rely heavily on Morse theory for 
the distance function p„. In the following sections we provide the framework and Morse 
theoretic ideas that are needed to prove each of the bounds. 


6. Upper bound - Critical points 

The main idea in proving the upper bound in Proposition 5.2 and Theorem 5.4 is the 
following. By Morse theory, changes to Hk occur only at critical points of index k and 
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k + 1. We will show that if A = ojdnr'^ — log n + k log log n + w(n) then there are no such 
critical points with value greater than r, and therefore the Hk{C{n, r)) = Hk{C{n, rmax))- 
On the other hand, with high probability we have that U{n,rjna,x) covers the torus, 
and therefore Hk{C{n,r)) = (cf. Corollary B.2). In this section we provide the 

definitions and statements needed to make this argument rigorous. 

The main object that will help us analyze the behavior of homology later is the 
following quantity - 

C'fe(r) ;= # {critical points c S of index k, with p„(c) < r} , (6.1) 

i.e. the number of critical points of index k with value in [0,r]. 


6.1. Expected value 

The following proposition provides the expected value of Ck{r), as a function of A. Note 
that as opposed to most of the results related to random geometric graphs and complexes, 
the following statement is true for every n and is not asymptotic. 


Proposition 6.1. For 1 < k < d, and for every r < rmax we have , 


E{C'/c(?')} = DkU 



where is a constant that depends on k and d and is given in (8.12), A is defined in 
(4.2), and rmax is defined in (3.1). 


An immediate corollary of Proposition 6.1 is related to the Euler characteristic of the 
Cech complex. The Euler characteristic is an integer valued topological invariant that can 
be viewed as a “summary” of the Betti numbers. There are various equivalent definitions 
for the Euler characteristic of a simplicial complex E. One of the definitions, via Betti 
numbers, is 

OO 

X(E) = ^(-l)'=/3fc(E). (6.2) 

k=0 

In other words, the Euler characteristic is a (signed) sum of the number of cycles of all 
possible dimensions. In the special case where E = C(n, r), defining 

X{r) ■■= x(C(n,r)), 

we can use Morse theory (see [34]) to show that 

d 

/c=0 


Thus, using Proposition 6.1 we have the following. 
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Corollary 6.2. If r < rmax then 


for some constants Aj that depend on j and d, defined in (8.14) below. 

Remarks: The expected Euler characteristic formula in Corollary 6.2 will not be used 
directly in proving the main results of this paper. Nevertheless, we included it here both 
for completeness, and for the fact that it provides valuable intuition which sometimes 
may even lead to concrete statements. At this point we want to mention the following 
observations. 

The terms Ajne~^A^ appearing in the sum above vanish at A — logn + j loglogn. 
This implies that a significant change in the topology of the complex should occur 
at these points. While this intuitive statement is somewhat vague, it does agree 
with and strengthens the main results of this paper. 

When A —)■ oo we have 

E{x(r)}«Ad_inA‘'-ie-^, (6.3) 

where Aa_i = and > 0. Thus, we can conclude that for very 

large A the Euler characteristic converges to x(T‘^) = 0, and its sign depends on 
the dimension of the torus. The intuition behind this phenomenon is that as we get 
closer to covering the torus two main things happen: 

• All the nontrivial cycles in the homology of have been formed, except for 
the cycle in the top dimensional homology. 

• The lack of coverage introduces many small voids (“bubbles”) generating non¬ 
trivial {d — l)-cycles. 

If K is the number of bubbles, the Euler characteristic would therefore be 
X(r) « x(T‘') - (-1)" + = {K + 

which explains why sign(x(r)) « (—1)“^“^ in this case. 

6.2. Vanishing thresholds for critical points 

Recall the definition of Ck{r) in (6.1), and define 

C'fc(r) := Ckirma^) - Ck{r), 

i.e. the number of critical points c with p„(c) G (r, rmax]- Since, with high probability we 
do not expect to have critical points with values higher than rmaxi the value of Ck{r) 


1 . 


2 . 


( d-l 

1 + X] 
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can be thought of as the “remaining” number of critical points. Using Proposition 6.1 
we have that 


E{C'fc(?')} = DkU 



„-A„ 



(6.4) 


where A = and A„iax = which implies that A ^ Amax (since r —>■ 0 and 

J’max is positive). If, in addition, we assume that A —>■ oo, the dominant part in (6.4) is 
Dkne~^A^~^, and therefore we have the following phase transition. 


Corollary 6.3. Suppose that A = logn + (fc — 1) log log n + w{n), then 

{ oo w{n) —)■ —oo, 

■u;(n) —)■ a S (— 00 , oo), 

0 w{n) —)■ oo. 

The above corollary suggests that A = logn + {k — 1) log log n is the threshold for 
the vanishing of Ckir), i.e. it is a point above which we do not have any more critical 
points of index fc. As mentioned in the beginning of this section, the vanishing of the 
critical points combined with Morse theory will help us prove the upper bound part of 
Proposition 5.2 and Theorem 5.4. 


7. Lower bound - ©-cycles 

To prove the lower bound part in Proposition 5.2 and Theorem 5.4 we will show that 
there exist nontrivial fc-cycles in Hk{C{n, r)) that do not belong to the torus T”* (i.e. they 
are mapped to trivial cycles in Hk{T‘^)). To prove that such cycles exist, we will be 
studying a special construction of fc-cycles that shows up when A —>■ oo, and is localized 
around critical points of index fc. 

Let V C and y G V. Recall that if y generates a critical point of the distance 
function then C{y) is the critical point itself, R{y) is the critical value (the distance), 
and A(y) is the embedded simplex spanned by y. For a fixed e S (0,1) define 

Ae{y) ■■= BR^y){c{y))\B:^^y^{c{y)), ( 7 . 1 ) 

where Br{x) is a closed ball of radius r centered at x, and B° {x) is an open ball, i.e. A^{y) 
is a closed annulus around C{y). Recall that by Morse theory, every index fc critical 
point either creates a new nontrivial fc-cycle or terminates an existing (fc — l)-cycle. The 
following Lemma provides a sufficient condition for the former to happen. While the 
statement is phrased for point sets in using the local isometry of the torus into , 
we can show that this is true for critical points c G T'^ with p^ (c) < rmax- 
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Lemma 7.1 (0-cycles). Let V C be a set of points, the distance function, and 
y CV be a set of k + 1 points generating a critical point of index k. Set 


cf = ^{y) 


infxgaA(y) Ik - 
2R{y) 


(7.2) 


where 9A(3^) is the union of the faces on the boundary of the simplex A(3^). If ^^(J^) C 
U{V,R{y)), then the critical point generated by y creates a new nontrivial cycle in 
Hk(C{V,R{y))). A cycle created this way will be called a “Q-cycle”. 


Remark: Note that for k = 1, i9A(3^) = y (the boundary of a 1-simplex are its end 
points). In this case we have (f{y) = 1 and then A^(y) = dB{y) (i.e. a (d — l)-sphere). 
The statement above is still true. 

The proof of this lemma (see Section 8) requires some more working knowledge of al¬ 
gebraic topology, than is required for the rest of this paper. Therefore, we wish to provide 
the following intuitive explanation. At a critical point of index k a new fc-dimensional 
simplex is added to the complex C(n, r), which is isolated in the sense that it is not a face 
of any higher dimensional simplex. We argue that adding this face necessarily generates 
a new fc-cycle. Since A^{y) is covered by U{V,R{y)), the (d — l)-sphere dB{y) is cov¬ 
ered as well, and so the new cycle is formed somewhere inside this sphere. The effect of 
adding the new A:-simplex can be viewed as “splitting” the void inside the sphere using 
a fc-dimensional plane. This action generates a nontrivial fc-cycle. For example, a line 
crossing a 2-sphere will generate a 1-cycle (or a loop), whereas adding a plane will split 
the sphere into two “chambers” forming a 2-cycle. Figure 4 provides visual intuition to 
this explanation. Note that in the case where d = 2 and fc = 1, we place an edge splitting 
the circle into two (creating a new 1-cycle), and thus forming a shape similar to the letter 
0, hence the name “0-cycles”. 

To prove that Hk{C{n, r)) ^ idfc(T‘^) we will show that /3fc(r) > /3fc(T‘*) by counting a 
subset of the cycles originating from 0-cycles. For a fixed e > 0, we define /3^(r) to be 
the number of 0-cycles generated by subsets y dVn such that 

(1) R(3^) e (r',r], (2) Br^,{Ciy))nrn=y, and (S) cf{y) > e, 

where 

r' = r{l — S), r" = r(l + V^), S = A~^. 

Condition (1) ensures that the 0-cycles are created at a critical radius in (r',r]. Con¬ 
dition (2) ensures that at radius r, the new fc-simplex added is still isolated, and therefore 
the new cycle it generates is nontrivial (this is justihed in the proof of the next Lemma 
7.2). Condition (3) is required in order to be at a safe distance from the singular case 
4>{y) = 0. In other words, PKr) counts 0-cycles that are created before r and terminated 
after r, as stated in the following lemma. 


Lemma 7.2. Let e > 0 be fixed, then 


Mr) > Mr)- 
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Next, we will evaluate the expected value of I3l{r). 

Lemma 7.3. There exists e > 0 such that i/A —^ oo, then 

E{/3|(r)}~nA'=-2e-A. 

This lemma provides the proof for the lower bound in Proposition 5.2. To prove the 
lower bound part of Theorem 5.4 we will need to prove the following lemma. 

Lemma 7.4. There exists e > 0 such that for any 1 < k < d — 1, if A = logn + (fc — 
2) log log n — w{n) and w(n) ^ log log log n then 

lim P {Pl{r) > 0) = 1. 

n—>-oo 



r < Riy) r = R(y) 

(b) 


Figure 4. A visual intuition for ©-cycles for d = 3. (a) A ©-cycle with k = 1: The set y consists of the 
two white dots placed on the sphere dB{y) of radius R{y). We drew the balls around the two points 
as semi-transparent and emphasized the disc at the equator to demonstrate the effect of the critical 
point. Recall that the assumptions in Lemma 7.4 imply that the sphere dB{y) is covered by 
Therefore, once the two balls touch (at r = R{y)) the gap between the balls and the covered sphere 
creates a closed loop, or a 1-cycle (marked in red). This effect is equivalent placing a 1-dimensional ball 
(a string) at the equator, (b) A ©-cycle with k = 2\ In this case the set y consists of three points. Once 
the balls around the points touch, they divide the void inside the sphere dB{y) into two, introducing a 
new 2-cycle. Topologically, this is equivalent to placing a 2-dimensional disc at the equator. 
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8. Proofs 

In this section we wish to present the detailed proofs for the statements in Sections 5- 
7. In Section 8.1 we prove the main results that appeared in Section 5, based on the 
lemmas that appeared later in Sections 6 and 7. Next, in Sections 8.2 and 8.3 we prove 
the lemmas that appeared in Sections 6 and 7, respectively. 

8.1. Proofs for Section 5 

Proof of Proposition 5.2. 

We need to show that aknh-^~'^e~^ < E{/3fc(r)} < /3fc(T‘*) + hknA^e~^. 

Upper bound: Denote /3/c(r) = /3fc(r)—/3/j(T'^), then we need to show that E{/3fe(r)} < 
hknK^e~^ for some bk > 0. This will be done using the expected number of critical points 
provided by Proposition 6.1. The main idea is that the number of critical points of index 
k + 1 occurring at a radius greater than r (i.e. C'fc+i(r)) should serve as an upper bound 
to $k{r), since all the extra cycles (the ones that do not belong to the torus) must be 
terminated at some point. 

Let E denote the event that fi(n,rniax) covers the torus. Then, 

E{/3fc(r)} = E{/3fc(r) | E}V (E) + E0k{r) \ E^V (E^). (8.1) 

If the event E occurs, then Pn(x) < r-max everywhere on the tours, and therefore using 
Morse theory (see Section 3.2) every extraneous A:-cycle in C{n,r) that does not belong 
to the torus will eventually be terminated by a critical point of index k + 1 with value in 
(j’Umax]- Therefore, we must have /3/c(r) < C'fc+i(r) = C'fc+i(rniax) - Ck+iir), and then 

E{/3fc(r) I E}¥ (E) < E{Ck+i{r) \ E}V (E) < E{C'fc+i(r)}. (8.2) 

Using Proposition 6.1 we have that 

( k ■ k \ 

’ (8.3) 

3=0 j=0 J 

where A = u!dnr‘^, and A^ax = oJd.nr'^^^, which implies that A <C Amax (since r —3 0 and 
J'max is positive). Since, in addition, we assume that A —> oo, the dominant part in (6.4) 
is Dk+ine~^A^. Putting it back into (8.2) we have 

E{/3fc(r) I U}P(£;) = 0{ne-^A^). 

To complete the proof we show that the second summand in (8.1) is o{ne~^A^). For any 
simplicial complex, the fc-th Betti number is bounded by the number of fc-dimensional 
faces (see, for example, [24]). Since the number of faces is bounded by (^j), we have 
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that 

E{/3fc(r) I E^}F (E^) < E | I ^ (E^) 

OO . ^ 

= E L“.P(IP-I = "'|S')P(E“) (8.4) 

7n—k-\-l ^ ^ 

= E (j^1^y(.E^\\'Pn\=m)¥i\rn\=m), 

m=k+l ^ ^ 

where we used Bayes’ Theorem. Since Vn is a homogeneous Poisson process with intensity 
n we have that P(|7^„| = m) = Jf , and also that given \Vn\ = m we can write Vn 
as a set of m i.i.d. random variables Xm = {-^i, ■ ■ • ,^m} uniformly distributed on the 
torus. Therefore, 

FiE^l |lP„|=m)=P(TVZ^(-^r„,w)). 

Let e = rniax/2, and let 5 be a e-net of i.e. for every x G T‘^ there is a point s G S 
such that p{s, x) < e. Since rmax is fixed, we can find such an e-net with |5| = M, where 
M is a constant that depends only on d. Note that if for every s G S there is Xi G Xm 
with p(s, Xi) < e, then for every a: G we have 

p{x, Xi) < p{x, s) + p{s, X^) < 2e = rinax, 

and therefore C ld{Xjn,rma,x)- Thus, if ^ ld{Xm,rmax), then there exists s G S 
with p^^ (s) > e, which yields 

P (T'" ^ U{Xm, w)) < ^ p {px^ (s) > e) = M(1 - ccde'^)’". 

s^S 

Putting it all back into (8.4) we have 

°° / rn \ p-n m n/r 

Since e > 0 is a fixed constant, and since A ^ n and A —)■ oo, we have that 

E{/3fc(r) I E‘^}F{E^) = o{ne-'^A^). 

Since /3fe(r) < /3fc(r) we also have 

E{/3fc(r) I £;^}P(£;^) = o{ne-^A^). 

To conclude, we have E{/3fc(r)} = 0{nA^e~^). That completes the proof for the upper 
bound. 

Lower bound: We need to show that E{/lfc(r)} > aknA^~^e~^ for some Ok > 0. 
From Lemma 7.3 we know that E{/3^(r)} ~ nA^~^e~^. From Lemma 7.2 we have that 
/3fc(r) > ld%{r). That completes the proof. □ 
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Next, we prove Corollary 5.3. 

Proof of Corollary 5.3. 

From Proposition 5.2 we have that E{/3fc(r)} < /3fc(T‘^) + bkne~^A^ and therefore if 
A = logn + fcloglogn + w(n) we have E{/3fe(r)} —>• /3fc(T‘^). In addition, we also have that 
^ Taking A = logn + (fc — 2) log log n — w{n) we then have that 

E{,5fc(r)} -)■ oo. □ 

Finally, we will prove Theorem 5.4. 

Proof of Theorem 5.f. 

Upper bound: From Corollary B.2 for every 8 G (rmax/2, Tmax) we have that C 
Z^(n, (5) with high probability, implying that Hk{C{n, 6)) = From Corollary 6.3 we 

know that if A = logn + A:loglogn + w(n) then both E{C'fc(r)} —)■ 0 and E{C'fe+i(r)} —)• 0. 
Using Markov’s inequality, this implies that with high probability both C'fc(r) = 0 and 
Ck+i{r) = 0. By Morse Theory, this implies that for 6 G (r,rmax) the fc-th homology of 
C(n, (5) does not change. Therefore, we conclude that with high probability Hk{C{n, r)) = 

Lower bound: If A = logn + {k — 2) log log n — w{n), and w{n) ^ log log log n then 
from Lemma 7.4 we have P{/3^(r) > 0) —>■ 1 for every 1 < k < d — 1. From Lemma 7.2 
we have that /3fc(r) > Pl{r). Therefore, with high probability we have /3fe(r) > 0 which 
completes the proof. 

□ 


8.2. Proofs for Section 6 

Proof of Proposition 6.1 . 

We want to show that E{C'fc(r)} = DkU ^1 — for r < rmax- First, note 

that the expected number of critical points with r = r^ax is zero, since this requires a 
subset of fc + 1 points to lie on a sphere of radius exactly rmax and this happens with 
probability zero. Thus, for the remainder of the proof we will assume that r < rmax- 
Let 7^ C be a finite set of points with its corresponding distance function 
and let c be a critical point of the distance function with p^ (c) < r. Following the 
discussion in Section 3.2 and Lemma 2.4, since r < rmax then c is generated by a subset 
y that satisfies the following: 

(1) C(3^) G A(3;), (2) B(37) n IP = 0, {3)R{y)<r, (8.5) 

where C{y), R{y), B{y) and A(3^) are defined in (2.4)-(2.7). For any x G and 

7^ C K’’*, we define the following indicator functions: 

fc(x) := 1 {C'(x) G A(x)} , 
hr{x.) := /i(x)l {i?(x) < r} , and 
Pr(x, V) := /ir(x)l {7?(x) n 7^ = 0} . 


( 8 . 6 ) 
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These functions can also be applied to points on the torus, provided that x is contained 
in a ball of radius less than rmax and using the isometry into If this is the case, we 
will define /i(x), hr.(x), and gr(x) the same way, otherwise we will set the function value 
to be zero. Using the above definition we can write 

Ck{r)= Y. 9r{y.Vn). 

ycv„ 

\y\=k+i 

Applying Palm theory to the mean value (see Theorem A.l) we have 

fc+l 

= (^TT)! ® ^ ’ 

where y' is a set of {k + 1) i.i.d. points on T'^, independent of Vn- Note that by the 
properties of the Poisson process Vn we have that 

E {1 {B{y) n = 0} I y'} = P (|U„ n B{y)\ = o | y) = e-’^Voi(B(y')), 

Since we require that R{y) < r < rmax then B{y) is isometric to a Euclidean ball, and 
Vol(B(y)) =WdA‘^(y). Therefore, 

Ai I J_ ^ 

= (8.7) 

[k + 1)! 

Note that the functions hr(x) and A(x) are translation invariant, and therefore we can 
write 

A) I 1 p 

(fc +1)! y(Td)fc 

where x' = (xi,... ,Xk}- Since the set (0,x') is contained in a ball of radius r, then 
necessarily we have that Xi S i?2r(0) for all 1 < i < fc and since r < r^ax we can 
isometrically embed this ball into Thus, from here on we can assume that we are 
integrating over points in and write 

Ai I X p 

(/c + 1)! J(Kd)fe 

To proceed, we use the following change of variables: 

Xi -)■ ryi, l<i<k, pi 

Denoting y = {yi,..., pk) £ we then have, 

dx' ^ 

y(0,x') hr(0,ry) = /ii(0,y), 
i?(0,x') A(0,ry) = ri?(0,y). 






24 


O. Bobrowski & S. Weinberger 


Therefore, 


E{Cfe(r)} 


j^k+l^dk 

(fc + l)! 





nA^ 

^1{k + iy. 


hi(0,y)e-^^'(°’y)dy. 


( 8 . 8 ) 


To further simplify this expression, we take the following change of variables. Instead of 
using the coordinates yi,... ,yk G we use the fact that the k+1 points - { 0 , yi,..., yk} 
lie on the unique (fc—l)-sphere ^(O, y) (defined in (2.3)). We can use spherical coordinates 
(r, 9) to identify the center of the sphere, where r G (0, oo), and 9 G (the {d— l)-unit 
sphere). Since 0 G S'(0,y) we have that i?(0,y) = r. Once we know where the center 
is, in order to uniquely identify the sphere S'(0,y), we need to choose the fc-dimensional 
subspace in which it resides. We represent this choice by the variable 7 G T^^. - an element 
in the fc-dimensional Grassmannian (the space of all fc-dimensional linear subspaces of 
M'^). Finally, after identifying the sphere 5(0, y), we need to choose k points on it, we do 
that by taking spherical coordinates - 9i,... ,9k G S^~^. Combined together we have the 
change of variables 

j/i-)■ r • (0 + 7 o 0*), l<i<k. (8.9) 

Defining J{t,9,j,0) := we can rewrite (8.8) as 


E{Cfc(r)} = 


nA^ 


^d(^ + l)-"'o Je,j, 


h{0,9 + ^o9)e J{T,9,^,6)d9d'yd6dT. 


Notice that (8.9) implies that J{t,9,^,9) = ^J(l, 0 , 7 , 6 ), and therefore 

nA^ ri , p 


E{Cfe(r)} = 


a;*(fc + 1 )! Jo 

/•i 


^dk-l^-Ar 


/ Oj-jjO 


h{0 ,0 + 70 9)J{1, 9, 7 , 9)d9d'yd9 


= DkuA^ 


^dk — 1 —At° 


dr. 


( 8 . 10 ) 


where Dk := fe-yO 9 + j o 9)J(1,9, 7 , 9)d9d'yd9 - a constant that depends 

only on k and d. Next, we make one final change of variables - Ar'^ —>■ t, which implies 
that dt = d ■ AT^~^dT and thus, 

E{C'fc(r)} = ^n f t^~^e~*dt. (8-11) 

“ Jo 

The latter integral is known as the lower incomplete gamma function and has a closed 
form expression which yields, 

-A 



E{C'fc(r)} = ^n{k - 1)! ( 1 - e 
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Finally, denotins 

p Dk{k-1)1 
d 

completes the proof. 

Remark: The conditions in Lemma 2.4 can also be phrased by saying that the points 
in y generate a critical point if their Voronoi cells have a common intersection with their 
dual Delaunay cell. The Voronoi and Delaunay tessellation have been extensively studied 
in the stochastic geometry literature. In particular, for Poisson-Delaunay complexes, it is 
known that circumspheres of the d dimensional cells (which correspond to the radius of 
index d critical points) have a gamma distribution [35] . This agrees with the form of the 
integral in (8.10), from which the incomplete gamma function in the Proposition comes 
from. 


( 8 . 12 ) 

□ 


Proof of Corollary 6.2. First note that since is a d-dimensional manifold and 
U{n,r) C T'^, then for k > d we have Hk{U{n,r)) = 0. For r < rmax, by the Nerve 
lemma 2.2 this is true for C{n,r) as well. Therefore, 


X{r) = ^(-l)'=/3fc(r). 


(8.13) 


k=0 


One consequence of Morse theory (cf. [34]) is that the sum of Betti numbers in (8.13) 
can be replaced by 

d 

k^O 

Using Proposition 6.1, and the fact that E{(7o(r)} = n for every r > 0 (all the minima 
are at distance 0), we have for every r < r^ax 


k-l 


E {x(r)} =n + J2i-1)'^ \Dkn\l- ^ 




k=l 


= n l + y](-l)'=Dfe 


J -0 

d-l 


.-A 




i=o k=j+l 


Denoting 


we have 


A = ^ E 


(8.14) 




d-l 


E{x(»')} = «(l-^o)+ne 

1=0 


(8.15) 
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Taking A —cx), from Corollary 5.3 we have that in this case E{/3fe(r)} —>■ /3/c(T'^) = 
Using (6.2) we then have 

hm E{x(r)} = x(T‘') = ^(-l)'=f^) =0. 

n—^oc) \ K / 

k^O ^ ^ 

On the other hand, as A —> oo (and n —> oo) the second term in (8.15) vanishes while 
the first keeps growing. In order for the limit to be zero, we therefore must have Aq = 1. 
That completes the proof. 

□ 


8.3. Proofs for Section 7 


Lemma 7.1 is central to the proof for the existence of cycles. However, its proof is more 
algebraic topological and in a different spirit than the rest of this paper. The details 
of this argument are not used elsewhere in the paper, and an intuitive explanation was 
presented in Section 7. 

In the following proof we will be using both simplicial homology (for the Cech complex 
C{V, r)) and singular homology (for the union of balls 7/(7", r)). An introduction to both 
can be found in [24, 36]. In both cases, we denote 

• Cfc - the chain group, i.e. all formal sums of simplexes (should not be confused with 
Ck{r) used in the paper), 

• dk ■ Ck ^ Ck-i - the boundary operator, 

• Zk - the cycles groups consisting of all chains j G Ck such that 9^7 = 0, 

• Bk - the group of boundary chains, i.e. all the chains 'j G Ck such that 7 = dk+ij 
for some 7 G Ck+i, 

• Hk - the A:-th homology group defined by Hk = ZkjBk- 
Finally, recall the definitions in (7.1) and (7.2) - 


and 


A,(y) := BR^y^[c(y))\B:^^y^[c{y)). 


<t> = 


infxGaA(y) llx- C'(3^)|| 

2R{y) 


Proof of Lemma 7. 1 . 

We want to show that for a critical point of index k generated by y, if Afj){y) C 
U{V,R{y)), then a fc-cycle is created. Consider the Cech complex Cr ■= C{V,r). Since 
R{y) is the critical value, then at r = R{y) a single new fc-simplex is added to the 
complex, call it A. Let r“ < r be a radius close enough to r such that C^- = Cr\^ and 
in particular such that 9A G Cr- ■ Recall also that for a critical point to exist we require 
that B(y) nP = 0, and therefore A is not a face of any (k + I)-simplex in Cr (by the 
construction of the Cech complex). 



Vanishing of Homology in Random Cech Complexes 


27 


The boundary dA is a cycle in Zk-i{Cr-)- Suppose that dA G Bk-i{Cr-), then 
dA = ^7 for some 7 G C'fe(Cr-). Taking 7 ' = A — 7 then d^' = 0, so 7 ' G Zk{Cr). 
However, since A is not a face of any (fc + l)-simplex in Cr, we must have 7 ' ^ Bk{Cr) 
and therefore 7 ' is a nontrivial fc-cycle in Hk(Cr). This cycle is not homologous to any 
cycle in Cj.- since it includes A. Thus, to complete the proof all we have to show is that 
dA G Bk-i{Cj.-). 

Recall that every A:-simplex in the Cech complex Cr is of the form a = [xq, ... , 2 ;^], 
where Xi GV and 

Brixo) n • • • n Br{Xk) 7 ^ 0 . 

Since we only consider k < d, we can embed every abstract fc-simplex a G back into 
a A:-simplex in and denote its embedding by 7r(cr). Setting Ur := U{'P, r) C there 
is a natural map from Ck{Cr) - the group of simplicial chains - to Ck{Ur) - the group 
of singular chains - that comes out from the Nerve Lemma 2.2 (using the embedding 
tt). Let Ik ■ Ck(Cr-) -G CkiUr-) be the natural map between the simplicial and singular 
chains. The induced map hk '■ Hk{Cr-) -G Hk(Ur-) is an isomorphism by the Nerve 
Lemma 2.2. If we can show that hk-i{dA) = 0 G Hk-i{Ur-) then that will imply that 
dA = 0 G Hk-i{Cr-), and we are done. 

Finally, note that since we assume = Arj,{y) C Ur, and since A 2 ^ C A^, we 
can choose r' to be close enough to r so that A2(j, C Ur'. Our choice of (j){y) implies 
that 7r(i9A) C H20, and therefore ik-iid/A) G C'fc_i(H20). Since A24, is a d-dimensional 
annulus, we have L7fe_i(^2<>) = 0, and therefore ik-i{dA) G Bk-i{A 2 ,j,) C Bk-iiUr-). 
Thus hk-i{dA) = 0, which completes the proof. □ 

Next, recall the definition of /3l{r) in Section 7, and define - 

huy) := {hr{y) - hr'{y))i wy) > e} , 
g^riy,v) := hiiy)i {Br"{c{y)) nVn = y}i {^,( 3 ;) C u{v,R{y ))}, 

where 

r' = r{l — S), r'' = r{l + V^), 6 = A~^, 

and hr is defined in ( 8 . 6 ). Then we have 

Pk{'r)= driy^'Pn). 
ycVr 


Proof of Lemma 1.2. 

We need to show that /3fc(r) > I3%{r). If gr{y,'Pn) = 1 then according to Lemma 7.1 we 
have that y creates a 0-cycle at radius R{y) G {r',r]. To complete the proof we will 
show that this ©-cycle created prior to r still exists at r and therefore contributes to 
Pk{r). 

In the proof of Lemma 7.1 we saw that when a 0-cycle is created a new fc-simplex A 
is added to the complex. This simplex is “isolated” in the sense that it is not a face of 
any (fc -|- l)-simplex. As long as A remains isolated, the cycle it creates cannot become 
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a boundary (a homologically trivial cycle). Thus, it remains to show that gliy^Vn) = 1 
implies that A is isolated at radius r. 

For the fc-simplex A generated by a set y to be isolated at radius r, we need to verify 
that nBr{y) := Br{Y), will not intersect with any of the other balls of radius 

r. We will show in Lemma 8.1 that the distance between points in r)Br{y) and C{y) 
is bounded by \/r‘^ — R?{y) < \/r‘^ — < rV^- Thus, by the triangle inequality, if 

the ball of radius r" = r(l + V^) around C{y) contains only the points in y, there 
is no other ball of radius r that intersects r\Br{y). Since g'r{y,'Pn) = 1 requires that 
Br"{C{y)) DVn =y 'we completed the proof. 

□ 


Lemma 8.1. For every x S nBr.(3^) we have \\x — C'(3^)|| < ^Jr"^ — E?{y). 

Proof. Let y = {ccq, ... ,Xk}, c = C{y), and R = R{y), so that ||c — Xi\\ = R for all 
i. Since c € A(3^) we can write c = where Oi > 0 and = 1. Let 

X S r\Br{y) then ||a; — a:i|| < r for every i. Now, 

\\x - x^\\^ = \\x - c + c - X.,\\^ = ||x - c\\^ + \\c - x^\\^ - 2 {x - c,c - x^), 

where (•, •) is the inner product in Thus, we conclude that for every i 

\\x — c\\^ < r"^ — R^ -\-2 {x — c,c — Xi). 

Summing over all the points Xi, we have 

k k k 

cti ||a: - c||^ < - i?^) + 2 Oi (a; - c, c - Xi). 

i—0 i—0 i—0 

Since Q^i = 1 and o^iXi = ewe have 

||a; — c||^ < — R^. 


□ 


Next, we prove Lemma 7.3. 


Proof of Lemma 7. 3. 

We need to show that E{/?^(r)} ^ We start by fixing e > 0 and evaluating 

the expectation of I3%{r)^ similarly to the computations in the proofs of Proposition 6.1. 


E{/3^(r)} 




(fc + l)! 

nk + l 


E{5^(y,y uiP„)} 


(A: + l)! 




(8.16) 
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where 


p,(x) := P {A^iy') C U{y U Vn, R{y)) I y = X, n S,yc(x)) = 0). 

We will first show that Pe(x) —1 uniformly for all x. Denoting P^, (•) := P (• | n (C'(x)) = 0), 
we observe that 

Pe(x) >PJA,(x) Cf^(P„,i?(x))). 

In the following we use the shorthand notation: 


R = R{x),A^ = A^{x),B = Br"{C{x)),p^ =Pe{^)- 


Let 5 be a (ei?/2)-net of Aj, i.e. for every cc S there exists s G S such that ||a: — s|| < 
eRl2. Note that there exists c > 0 such that 


Vol(A) uJdR‘^(l-e‘^) 

\S\ < C -,—= c , V ^ = cu;d (2 efn - e‘^) =: Ci. 

' ' “ {eR/2Y (eB/2)^ ay / j y j 

By the triangle inequality, if for every s G S we have Vn H Bii(^i_^/ 2 ){s) ^ 0, then 
A^ C U{Vn, R)- Therefore, 


Pe > PJVs G 5 : Bn n Bi^(i_,/2)(s) ^ 0). 

Note that under P^ we have already B VVn = ^- Thus, for any s G S 
V.i'Pn n Bdiil-e/ 2 ){s) = 0) = 

Recall that r" is the radius of B, and notice that since R > r', we have that r” /R < 
r"jr' -G 1. In particular, for n large enough we have r" < B(1 + e/4) (for all x). Since 
s G A^, then at least a fixed proportion of the volume of Bfl(i-e/ 2 )(s) must be outside 
B such that 

Vol(Bfl(i_e/ 2 )(s)\B) > C2UJdR‘^, 
for some C 2 G (0,1). Therefore, 

P,(3s GS:r„n Bfl(i_,/ 2 )(s) = 0) < ^P,(P„ n Bn(i-e/ 2 )(s) = 0) < 

sGS 


For every x such that h/(x) ^ 0 we have i?(x) > r', and therefore 1 — < 

Pe(x) < 1, Since uJdn{r'Y —t oo, we have Pe(x) —1 uniformly for every x. Going back to 
(8.16) we therefore have 

„A ;+1 r 

(fc + 1)! J (jd-^k+i 

Using the same change of variables we had in the proof of Proposition 6.1, similarly 
to (8.10) we have 

E{/3fe(r)} « 

Jl-S 
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where 

Dl := fc ^ [ h^O,0 + 'Yoe)J{l,9,j,e)d0d'rde, 

uj^{k +1)1 

and := h{y)l. {<p{y) > e}. Recall that the points in the last integral are arranged on 
a sphere of radius 1 centered at 0. Thus, (7.2) turns into 


^( 0 , 0 + j o 6) 


inf llx — 011 . 

xGdA{0,9-\-'yo6) 


This function is continuous in 0, 7 ,0, its minimum value is zero and it is achieved 
when the center of the sphere (0) lies on one of the faces of the simplex A(O,0 + 70 
9). This set has measure zero, and therefore there exists €k > 0 such that the set 
{ 0 , 7,0 : (j){O,0 + 7 o 0) > ek} has a positive measure. Taking e = €k we have that Df. > 0 
and since T‘^^~^dT = S + o(S) we have 




(8.17) 


Recall that 6 = A then 


uMr'Y = A(r"/r)‘^ = A (l + = A (l + dV^ + o ) 

Therefore, 

E{/3^(r)}«7?^.nA'=-2e-^-0(i). 

So 

E{/3|(r)}~nA'=-2e-A. 


That completes the proof. 

We will end this section with the proof for Lemma 7.4. 


A + 0(1). 


□ 


Proof of Lemma 

From Lemma 7.3 we have that if A = log n + (fc — 2) loglogn —?ii(n) then E {/3|(r)} —>■ 00 . 
If we show in addition that Var (/3^(r)) <C E {ld%{r)Y, then using Chebyshev’s inequality 
we have P(/3^(r) > 0) —>■ 1 (cf. Chapter 4 in [3]). We therefore need to find an upper 
bound on the variance. Recall that 

Var(^^(r)) = E{(^^(r))2} - E{/3^(r)}2. 

Since PI = J2ycv„ we can write {Pl(r))‘^ as 

mr))^= E E 9r{yu'Pn)9Uy2,Vn) 

yicVr, y2CVr. 
fc + 1 

= E E 9r{yi,'Pn)gt{y2,'Pn), 

j=o \yiny2\=j 
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where 3 ^i, 3^2 iterate over all subsets of Vn with {k + 1) points. 

Denoting the inner sum by Ij, observe that Ik+i = Pkir), and therefore we have 

k 

Var(/3^(r)) = E{/3^(r)} + ^E{/,} + (E{/o} - E{/l^(r)}2). 

To complete the proof we therefore need to show that E{/j} /E{/3^(r)}^ —>■ 0 for every 
1 < i as well as (E {/g} — E{/3|,(r)}^)/E{/3^(r)}^ —>■ 0. 

For 1 < j < k, using Palm theory (see Corollary A. 2) we have that 

2/-_|_2_ j 

= j!((r+i-j)!)2 ^ {/(X, y u rn)g^Ay2,y 

where X, 3^2 are sets of fc + 1 random variables, \y[ C 3 ^ 2 ! = h such that 3^' = 3’( U 3^2 is 
a set of (2A: + 2 —j) i.i.d. random variables, uniformly distributed in and independent 
of Vn ■ Since r < rmax we have, 

2^—j—2 J p 

-w/ h^(xi)h^(x2)e-"^("-’^^)dx, (8.18) 

j!((fc + 1 - j)!)^ 7(Td)2fc + 2-i 


where 


x= (a;o,---,a; 2 fe+i-j) G ^ 

xi = {xo ,... ,Xfc) G (T‘*)''+\ 

X 2 — (^ 0 ; ■ • • ; — 15 ^fc+1 J • ■ ■ J ^2k + l—j ) G (T ) , 

F(xi,X 2) = Vol(B(xi) U B{k2)). 

Note that (8.18) is an inequality, since we dropped the term verifying the coverage of the 
Aj-s, as well as the term verifying that xi C i?(x 2 ) = 0 and vice versa (both terms are 
bounded by 1). Notice that 


C(xi,X2) =uJdR'^{xi) +uJdR‘^{x2) - Vol(B(xi)nB(x2)). 


Let 14 (A) denote the volume of the intersection of two balls of radius r such that their 
centers are at distance A from each other. Since i?(xi) < r we have 

Vol(B(xi) n i?(x 2 )) < 14(||C'(xi) - C(X2)||). 

In addition, we know that i?(xi) > (1 — 6)r, and therefore, 

^2k-\-2—j 


E{/,}< 


n 


j\{{k+ i-j)\Y 

^2k+2—j 


= -A( 2 ( 1 -< 5 )'‘- 1 ) 


„-A 


j!((fc + 1 — j)!)2 Jfjd-^2k+i-j 


/ h^(xi)/i^(x2)e-^+^’-(ll‘^(’^i)-^(^^)ll)d> 
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since S = A Since the sets Xi, X 2 have at least one point in common (xq), we can apply 
a similar change of variables to the one used in the proof of Proposition 6.1 and have 


E{/,}< 


^2k+2-j ^d{2k+l-j) 

j!((/c +1 - j)!)2 


,-A 


/(Kd)2fc + 1 




where 

y=(yl,...,2/2fe+l-,)G(M‘'^+'-^ 

yi = iyi,---,yk) e 

y2 = (yi, ■ • ■ ,yi-i,yfe+i, • ■ ■ ,y2fc+i-i) g 


Let S' := {y : /if(0,yi) = hl{0,y2) = 1}, then S is bounded and we can write 


Js 

for some constant ci > 0. Next, fix a > 0, and define 


51 = {yGS:||C'(0,yi)-C(0,y2)||<a}, 

5 2 = S\Si = {y G S : ||C'(0,yi) - C(0,y2)|| > a} . 


We will bound the integral on both sets separately. On the set Si, note first that 
^_A+nrWii\\C{ 0 ,yi)-Ci 0 ,y 2 )\\) < ^nd therefore 


/ g-A+rirWi(||C(0.yi)-C(0,y2)||)^y 

Jsi 

<[ hl{0,yi) f /i"i( 0 ,y 2 )l {||C(0,yi) - C'( 0 ,y 2 )|| < a}dy. 

Jyi Jy2\yi 

For a fixed yi, if 11(7(0,yi) — (7(0,y2)|| < a, and i?(0,y2) G (1 — i5,1) then for every 
y G y 2 \yi we have 

l-(5-a< ||y-(7(0,yi)|| < 1 + a. 

Thus, the (fc + 1 — j) points in y 2 \yi lie in a thin annulus, and we have 



^)(0,y2)l {110(0,yi) - (7(0,y2)|| <a}d(y2\yi) < {ujd{{l + af 




« (a;d(2da + 


where we used the approximation (1 + x)‘^ « 1 + dx (which will be justified later when 
we take a —>■ 0). Since i?(0,yi) G (1 — <5,1) we have that df(0,yi) = 0{6) = 0{A~‘^) 
and therefore. 


f g-A+rirWi(||C(0,yi)-C(0.y2)||)^y < C2A~^{2da + dS)’^+^-^ 
Jsi 


( 8 . 20 ) 
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for some constant C 2 > 0 . 

On 5*2 we use the fact that (see Appendix C): 

Vi(A) = ujd — ujd-iA + o(A), 

and together with the fact Vi is decreasing in A, we have 

f g-A+nr‘*yi(||C(0,yi)-C(0,y2)||)^y < < Vol(S'2)e 

JS2 

for some C 3 > 0. Combining (8.19), (8.20) and (8.22), we have that 

E{Jj} < C4nA2'=+i-^e-^(A-2(2a + (5)'=+i-^' +e-"^^“), 
for some constant C 4 > 0. From Lemma 7.3 we have that 

E{/3^(r)} ~ 

and therefore 

For A = log n + (fc — 2) log log n — w{n) we have 

EU 4 } 


( 8 . 21 ) 

— csAck 

5 

( 8 . 22 ) 


< 


((logn)'=+i-^( 2 a + + (logn) 


for some C 5 > 0. Taking 

fc + 3 — j log log n 

a = -^-, 

C 3 log n 

and assuming that w{n) ^ log log log u we have 

E{/,} 

E{/3^(r)}2 “ ■ 

To complete the proof, it remains to show that 

E{Jo}-E{/3^(r)}2 


fc+3-i -C 3 log(n)Q 


)■ 


E{/3^(r)}2 


^ 0 . 


Using Palm theory, we have that 


E{Jo} = 


^2k+2 


{{k + i)\f 


E{ff^(X,y U Vn)gt-{y'2,y U Vn)} , 


where y, y are two disjoint sets of fc + 1 i.i.d. random variables, uniformly distributed in 
independent of Vm and 3^' = y U y. On the other hand, we can rewrite E{/l^(r)}^ 
as 


E{/ 3 ^.(r)r = 


, 2fc+2 


{{k + v)\y 


rE{g^(y,y uu„)5^(y,yuu;)}, 
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where the same as above, and is an independent copy of Vn- Next, define 

Ag := g^riy'i,y' u Vn)9i{y2.y u vy - 9i{y[,y[ u Vn)9i{y2.y2 u K) 

then 

„ 2fe+2 

E{Io} - Emir)r = ^^^—^E{Ag} . 

We can split Ag into two terms - 

Ai :=Ag-l {B{y[) n = 0} A 2 := Ag ■ 1 {B{y[) n ^ 0} . 

First, we show that E{Ai} = 0. Observe that 

Ai = K{y[)K{y^)\{B{y[) n B{y^) = 0} 

X (i {Vn n B{y[) = 0} 1 [Vn n Biy'^) = 0} 

- 1 {ip„ n B(X) = 0} 1 {iP; n = 0}). 

If Ai 7 ^ 0, then B{y[) and Bly'^) must be disjoint. Therefore, given y{ and y!^^ the set 
Vn n 5 ( 3 ^ 2 ) is independent of the set Vn H B{y{) (by the spatial independence of the 
Poisson process), and has the same distribution as VyB{y 2 ). Thus, E {Ai | 3^2} = Oj 

for any 3^i,3^2- This implies that E{Ai} = 0. For A 2 , notice that 

A 2 < gl{y[,y[ U Vn)gl{y 2 y 2 u ^n)i {S(X) n Biy'^) + 0}, 


and then 


E{A2} < [ /ir(xi)/i^(x2)e {B(xi) n i?(x2) 7^ 0} dxidx2, 

J(Jd'f2k + 2 

where xi = (xq, . •. ,Xfe),X 2 = (x^+i,... ,a; 2 fe+i). Note that here the sets xi,X 2 do not 
share any point in common. However, since the balls intersect, the points have to be very 
close to each other, and in particular to Xq. Thus, we can apply similar arguments to the 
ones we used for Ij that to show that 


lim 

n—^oo 


E{Jo}-E{^^(r)p 

E{/3^(r)}2 


= 0 . 


That completes the proof. 


□ 


9. Riemannian Manifolds 

The arguments in this paper should extend to all compact Riemannian manifolds. In 
this section we wish to sketch a heuristic argument for this. We hope to give a detailed 
treatment (with more effective results) elsewhere. 
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For flat tori, we niade use of the fact that for r sufficiently small, the r-balls are all 
(1) embedded (i.e. are topological balls) (2) flat (isometric to Euclidean balls), and when 
a sufficiently large number of points are sampled uniformly - (3) all the critical points 
of the distance function are associated with points within distance r of one another. We 
have also used the fact that (4) intersections of balls at small scales are either empty or 
contractible. 

For general compact Riemannian manifolds (1) and (3) are both true and straightfor¬ 
ward. (4) is also true when r is smaller than the convexity radius. However, (2) is not 
true in general, because the curvature might be nonzero. 

However, the problems that we are considering here are asymptotic and scale invariant. 
This formally implies that flatness is irrelevant; curvature can always be assumed to be 
arbitrarily small. 

More explicitly, scale r on the Riemannian manifold (A4,g) is the same as scale 2r on 
the manifold (A4, 2g) whose curvature tensor is multiplied by 1/4. Moreover, the covariant 
derivatives of the curvature tensor are multiplied by even smaller negative powers of 2. 
As a result, by working with a small enough r, we can arrange (in our original manifold) 
for the set of r-balls to be a precompact set in any fixed C^-topology (given by a bound 
on the {k + l)-st derivatives of the metric) in a small neighborhood of the flat gij = 6ij 
metric. A continuity argument then gives that for n large enough (so that all the critical 
points are at distance smaller than r) all of the integrals involved in the computations of 
expected values and variances are as close as desired to the flat values. The arguments 
of the paper then directly apply. 


10. Conclusion 

In this paper we studied the phase transition describing the vanishing of homology in 
random Cech complexes constructed over the d-dimensional torus. We saw that, as op¬ 
posed to other types of random complexes, for the Cech complex there exist only two 
sharp phase transition - one for Hq and one for the rest of the homology groups Hk, 
k>l. 

A possible explanation for the difference between connectivity and higher-order ho¬ 
mology is the following. To show that the Cech complex is connected, it is enough to 
show that there are no isolated points (See [40]). A point in the Cech complex is isolated 
if and only if the ball of radius 2r around it is empty (i.e. containing no other points 
from Vn)- On the other hand, as we have seen in this paper, the vanishing of the rest 
of the homology groups is related to critical points, and a critical point with value r 
requires that the ball of radius r around it is empty. In other words, connectivity is re¬ 
lated to empty balls of radius 2r whereas higher homology is related to empty balls of 
radius r. Heuristically, this should imply that the connectivity radius should be half the 
radius needed for the vanishing of the higher homology. In fact the difference between 
Aq = 2“‘^logn and A^ « logn describes exactly that (recall that A = Udur^). 

In addition to the sharp phase transition at A = log n we proved that there is a lower 
order scale such that the different homology groups vanish in an ascending order. The 
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results in this paper show that for Hk the vanishing threshold is in the range [logn + 
(k — 2) log log, logn + A: log logn], and it remains a future work to determine the exact 
vanishing point within this range. 
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Appendix 

Appendix A: Palm Theory for Poisson Processes 

The following theorem will be very useful when computing expectations related to Poisson 
processes. 

Theorem A.l (Palm theory for Poisson processes). Let {X,p) he a metric space, f : 
X —> K 6e a probability density on X, and let Vn be a Poisson process on X with intensity 
= nf. Let h{y, X) be a measurable function defined for all finite subsets y <Z X C. 
with |y| = k. Then 



yaVr, 

where y' is a set of k i.i.d. points in X with density f, independent ofVn- 

For a proof of Theorem A.l, see for example [40]. We shall also need the following 
corollary, which treats second moments: 

Corollary A.2. With the notation above, assuming l^ij = |y 2 | = k. 





13 ^ 10 ^ 21=1 

where y' = y'lC yl^ is a set of 2k — j i.i.d. points in X with density f, independent of 
Vn, and \y[ny!,\=j. 

For a proof of this corollary, see for example [9]. 

Appendix B: Coverage threshold 

The work in [20] studied the number randomly placed balls of radius r needed in order 
to cover each point in a compact Riemannian manifold m times. Translating the results 
in [20] for m = 1 using the notation in this paper, we have the following. Let be the 
number of balls of radius r required to cover a manifold M (assuming that the centers 
are uniformly distributed in Ai). Let a = uJdc'^ = n~^K, and define 


Xj. = aNn— \-Og{a ^) — dloglog(a ^). 

Theorem B.l (Theorem 1.1 in [20]). There exists ri > 0 and C > 0 such that 

F(W > a;) < Va;>0,r<ri, 

P (Xr < a;) < Ce"^, Va;<0,r<ri. 


(B.l) 

(B.2) 
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We wish to use this result to prove the following phase transition. 

Corollary B.2. Let M. be a compact Riemannian manifold, let Vn be a homogeneous 
Poisson process on A4 with intensity n, and let w(n) —>■ oo. 

1. //A = logn + (d — 1) loglogn + r(;(n) then 

M CU{Vn,r) w.h.p. 

2. If A = log n+ {d — 1) log log n — w{n) then 

M U{Vn,r) w.h.p. 

Proof. Recall our definition of the Poisson process in Section 4.1 as Vn = {Xi, • ■ •, ATat} 
where N ^ Poisson (n). For any n > 1 (not necessarily integer) define Xn = {-Ai,..., X^n^ 
and let A = ujdhr’^. 

1. Note that P {Xr > a;) = P (Ai U[Xn,r)), where x = an—log(a“^)—dloglog(Q;“^). 
Now, taking A = logn + (d — 1 ) loglogn + w{h), since a = ujdr‘^ = n~^A we have 

a; = (d — 1) loglogn + log A — dlog(logn — log A) + w{n). 


Note that 

lim (log A — loglogn) = lim (loglog n — log(logn — log A)) = 0. 

n—>-oo \ / n—>-oo \ / 

Therefore, we have x = w{n) + o(l), and using (B.l) we have 

P (A4 ^ U{Xn, r)) < (B.3) 

Moving back to the Poisson process, let q G (1/2,1) and set n = n — n^. Then 


'(MifIl(Vn,r))=P(M^Il(Vn,r) | |iP„| > A) P (|iP„| > n) 
+ P(MifIl(Vn,r) I IVnI <fl)P(IVnl <ii) 
<P(M^Il(Xn,r))+P(lVnl < A). 


(B.4) 


Using Chebyshev’s inequality we can show that P(|7^„| < n) —>• 0. Now, take A = 
logn + (d — 1) loglogn + w(n), and notice that A := ^A = A — Since 

lim (logn — logn) = 0 , and lim (loglogn — loglogn) = 0 , 

n—>-oo n—>-oo 

we have that 

A = logn + (d — 1) loglogn + w(n) + o(l). 

Combining (B.4) and (B.3) we then have 

P(Mif U{Vn, r)) < + o(l). 


Taking w{n) —>■ oo completes the proof. 


Vanishing of Homology in Random Cech Complexes 


39 


2. Notice that P {Xr < a;) = P (A 4 C U{Xn, r)) where x = ah—log(a“^)—dloglog(a“^). 
Similarly to the above, taking A = logh + (d — 1) log log h — w{h) and using (B.2) 
we have that 

CUiXf^r)) < 

Next, set h = n + for some q S (1/2,1), then 

P(M cZ^(lP„,r)) =P(M cZ^(iP„,r) | |iP„| > h) P (|iP„| > h) 

+ F{McU{Vn,r)\ \rn\<fi)F{\Vn\<fi) 

<F {\Vn\ > h) +F {M CU{Xf„r)). 

Continuing the same way as above yields 

P (Ad c UiVn, r)) < + o(l). 


and taking w{n) —>■ oo completes the proof. 


□ 


Appendix C: Intersection of Balls 

In this section we provide a brief explanation to the approximation of the volume Vi(A) 
that appeared in (8.21). 

The intersection of two balls is a region bounded by the union of two spherical caps. 
In [30] an explicit formula is shown for such volumes, from which we obtain 

^cos“^(A/2) 

Vi(A)=2a;d_i / sin‘^(6*)dd. 

Jo 

Thus, 

Vi'(A) = 2ujd-i — . ^ — sin‘^(cos~^(A/2)). 

’ 2^1- (A/2)2 ^ ^ ^ 

At A = 0 the intersection is a whole ball, and therefore Ti(A) = 0, and from the last 
formula we have V/(0) = —utd-i- Using taylor expansion we therefore have 


Ui(A) — ojd — u!d-iX + o(A). 
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